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Abstract 

The Dirac equation, with position-dependent mass, is solved approximately for the generalized 
Hulthen potential with any spin-orbit quantum number k. Solutions are obtained by using an 
appropriate coordinate transformation, reducing the effective mass Dirac equation to a Schrodinger- 
like differential equation. The Nikiforov-Uvarov method is used in the calculations to obtain energy 
eigenvalues and the corresponding wave functions. Numerical results are compared with those given 
in the literature. Analytical results are also obtained for the case of constant mass and the results 
are in good agreement with the literature. 
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I. INTRODUCTION 



In recent years, following work done by von Roos, and Levy-Leblond [1, 2] on relativistic 
and non-relativistic motion, the position-dependent mass (PDM) formalism has recieved 
attention in quantum physics. The solution of the Schrodinger equation with effective mass 
is a useful starting-point for the investigation of some physical systems, such as the electric 
properties of quantum wells and quantum dots [3], impurities in crystals [4-6], and semi- 
conductor heterostructures [7]. In general, for both relativistic and non-relativistic cases, 
the energy spectra and corresponding eigenf unctions have been studied by many authors 
using different methods and potentials: deformed algebras applied to a Coulomb problem 
[8], analysis within supersymmetric quantum mechanics [9-11], point canonical transforma- 
tion to study different classes of potentials [12], non-relativistic Green's functions applied to 
the harmonic oscillator [13], the Dirac equation in the Coulomb field [14], hyperbolic-type 
potentials [15], the Morse potential [16], and the Scarf II potential [17]. 

Another interesting area that has received a lot of attention is the investigation of so- 
lutions to the Dirac equation when the Dirac Hamiltonian has either spin or pseudospin 
symmetry [18]. The Dirac Hamiltonian of a particle with mass ji moving in scalar, V s (r), 
and vector, V v (r), potentials is invariant for two cases: the sum or the difference of the 
potentials is nearly zero. Important features of these symmetries provide an understanding 
of the structure of deformed nuclei, including deformation and super deformation [19, 20], 
and enable construction of an effective shell-model coupling scheme [21-28]. 

Recently, pseudospin and spin symmetry have been studied for various potentials, such 
as the Morse potential [29-31], the Woods-Saxon potential [32], the Coulomb potential [33], 
and the harmonic potential [34-36]. Some authors have also solved the Dirac equation in the 
context of pseudospin symmetry under the effect of the Eckart potential [37, 38], and the 
Poschl- Teller potential [39, 40], for spin-orbit quantum number k — 1 and/or any k- value. 
In Ref. [41], the bound states of both the Klein-Gordon and Dirac equations for the Hulthen 
potential are studied by writing these equations as a Riemann-type equation. In Ref. [42], 
the energy spectra of the Dirac equation for equal scalar and vector parts of the Hulthen 
potential is studied using a perturbative approach. In addition, the Dirac equation is solved 
for the Hulthen potential, for both spin and pseudospin cases, using the asymptotic iteration 
method [43]. 
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In the present work, we solve the effective mass Dirac equation for the generalized Hulthen 
potential within the framework of an approximation to the k(k — l)/r 2 term. We give the 
energy eigenvalue equation and the corresponding eigenfunctions for any spin-orbit quantum 
number k. We also obtain the energy eigenvalue equation in the case of constant mass. We 
give separately the Dirac solutions for the cases where spin and pseudospin symmetries are 
taken into account. We apply the parametric generalization of the Nikiforov-Uvarov (NU) 
method to obtain the energy eigenvalue equations and the eigenfunctions of the generalized 
Hulthen potential. The NU method is a powerful method for solving second order differential 
equations [44], like methods based on Lie algebras [45]. In Ref. [46], this approach was 
generalized to a theory based on spectrum generating algebras including relativistic effects. 
The NU method can be used, for both non-relativistic and relativistic CclSGS, clS el tool to 
find the energy spectrum and corresponding wave functions. 



II. DIRAC EQUATION 

The Dirac equation for a spin-| particle with mass ji moving in scalar V s (r), and vector 
V v (r) potentials is written as 



[a.P + (3[fx + V s (r)\ + V v (r) - E]^(r) = 0, (1) 

where E is the relativistic energy of the particle, P is three-momentum, and a and (3 are 
4x4 Dirac matrices [47] defined, respectively as, 



•-(::)■ •-(::) 

Here, a is a three-vector whose components are Pauli matrices, and / denotes the 2x2 
unit matrix. For a spherically symmetric potential, the total angular momentum J, and 
spin-orbit operator K = — /3(o\L + 1) commute with the Dirac Hamiltonian, where L is the 
orbital angular momentum operator. The eigenvalues of the operator K are k = ±(j + 1/2); 
k = -(j + 1/2) < stands for the aligned spin j = £ + 1/2, and k — (j + 1/2) > 
stands for the unaligned spin j = £ — 1/2. On the other hand, the total angular quantum 
number can be written in terms of the pseudo-orbital angular momentum £ = £ + 1 and the 



pseudospin angular momentum s = 1/2 as j = £ + s. The spin-orbit quantum number k is 
related to the orbital angular quantum number by the expressions k(k + 1) = £{£ + 1), and 
k(k — 1) = £{£ + 1) for a given k — ±1, ±2, . . .. The spherically symmetric Dirac spinor can 
be written in terms of upper and lower components, 



.«<r)=i( F ~ W1 V'-n (3) 

where Yj m (9, 0), and Y? m (0, 0) are spherical harmonics, and F nK (r)/r, and G nK (r)/r are the 
radial parts of the upper and lower components. Substituting Eq. (3) into Eq. (1) one can 
write the Dirac equation as a set of two coupled differential equations in terms of F nK (r) 
and G nK if) 



( t ~ ~r ) GnK (r) = [/i " EnK + V+ ^ F ™ ^ ' ® 
where V_(r) = V v (r) — V s (r), and V + (r) = V v (r) + V s (r). Using Eq. (4) for the upper 
component, and substituting into Eq. (5), we obtain generalized differential equations for 
the position-dependent mass case, 




k(k+ 1) 



- [/i + E nK - V-(r)] [pi - E nK + V + (r)} 



I dfj,(r) dV- (r) ' 

^ dr dr 



)(£ + 



pL + E nK -V_(r) 



\F nK (r)=0, (6) 



y dr 2 r 2 

( djMr) , dV+(r) \ / _d_ _ k\\ 

v dr i ^ dr rj_ \ fi ( \ n (n\ 

»-E nK + V + (r) ; G -^" U - ^ 
Here, the energy eigenvalues depend on the quantum numbers n, and k, and also on the 
quantum number £ according to the relation k(k — 1) = £(£ + 1). To solve these equations, 
we approximate the centrifugal term. Thus, the energy spectra and the corresponding 
eigenfunctions can be obtained analytically by using the NU method. 



III. PARAMETRIC FORMULATION OF NIKIFOROV-UVAROV METHOD 

We briefly give the mathematical background required for the parametric NU method 
[48], where the general form of the Schrodinger-like equation for any potential is written as 



Ms) = 0. (8) 



ds 2 s(l — 0:3s) ds [s(l — «3<s)] 2 
Comparing Eq. (8) with the general form of the basic equation of the method, we obtain 

f(s)=a 1 -a 2 s ; a(s) = s(l - a 3 s) ; a(s) = -£i« 2 + £ 2 s - £3 ■ (9) 



Substituting Eq. (9) into the polynomial 7r(s) = — 2 — ± J ( £ ~2 Z ) — o" + ka, we obtain 



7r(s) = o 4 + 0:5s ± \J (a~6 — ka 3 )s 2 + (o 7 + fc)s + o 8 , (10) 

where o 4 = \ (1 — ai), a 5 = \ (a 2 — 2a 3 ), o 6 = a!? + £1, o 7 = 20:40:5 — £ 2 , and o 8 = a\ + £3. 

In the NU-method, the function under the square root in Eq. (10) must be the square of 
a polynomial [48]. This condition gives the roots of the parameter k, 



h,2 = -("7 + 2o 3 o 8 ) ± 2^0309 , (11) 
with o 9 = 030:7 + o|o 8 + olq. For k = —(07 + 2a 3 o 8 ) — 2 v /o 8 og, n(s) becomes 

7r(s) = o 4 + a 5 s - [( V / «9 + Oi 3 y/o^ )s - y/o^} , (12) 

and also 

r(s) = «i + 2o 4 - («2 - 2o 5 )s - 2 [(v^ + <y 3y /a^)s - ^/o^} . (13) 

To satisfy the condition that the derivative of the function r(s) should be negative, we 
impose 
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r '( s ) = -(«2 - 2a 5 ) - 2( v /a^ + a 3 yo^) 
= -2a 3 - 2(^/a^ + a 3y /a^) < 0. 

The energy eigenvalue equation is written as [48] 



(14) 



a 2 n - (2n + l)a 5 + (2n + l)(y/ag + a 3y /a^) + n(n - l)a 3 

+ a-j + 2a 3 a 8 + 2* s /a%a$ = 0. (15) 

and the equality (x'(s)p(s) + <r(s)p'(s) = r(s)p(s) gives 

p(s) = s^-^l - q^)^"* 10 " 1 • (16) 
This equation together with y n (s) ~ ^j£^r [c n (s) p(s)] gives 

^( S ) = Pi ai0 ^-^ 10 " ) (l-2a 3S ), (17) 

where a w — a± + 2a 4 + 2^/as, an = a 2 — 2a 5 + 2(y / og + a 3 ^/a^) and P^ a '^(l — 2a 3 s) are 
the Jacobi polynomials. The equality = |^ [48] gives 

0( s ) = s ««(l-a 3S )- aia -^, (18) 
and the general solution ^(s) = (f)(s)y(s) becomes 

^ (s) = s « u( i _ a3 S )- ai2 -^pi ai ° _1 '^^ 10_1) (l - 2a 3 s) , (19) 
where oiyi = «4 + and 0:13 = 0:5 — (y^aig + a 3v /o8) [48]. 

IV. BOUND-STATE SOLUTIONS 

In order to solve the Dirac equation, including the term proportional to 1/r 2 , we set the 
vector and scalar potentials as the generalized Hulthen potential [31, 49] of the forms 
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W = V ° ~pi^—\ '^o>0, (20) 

-2ar 

V '( r ) = ~^o e _ 2ar _ 1 ,So>0, (21) 

where ct is the screening parameter (positive), the constant parameter S denotes the scalar 
and Vq denotes the vector part of the potential, respectively. Using the last two equations, 
we obtain 

-2ctr 

V-(r) = (Vo + S,)-^—^, (22) 

V + (r) = (Vo-So)^^- (23) 

where the "effective" potential V_(r) is a repulsive one in the relativistic region [18, 26-28, 
35, 50, 51]. It is well known that the Hulthen potential gives bound states in the non- 
relativistic region only if the potential parameters satisfy the condition that 5 2 < Vq < 45 2 , 
where 5 2 = 4a 2 , in the absence of the scalar part [41]. 

We use the equality dfx(r)/dr = —dV + (r)/dr to eliminate the last term in Eq. (7), and 
obtain the mass function which can be written as 

M0=*-3^. (24) 

The mass function has the same form as the Hulthen potential, where /xo denotes the 
integral constant, and /xi = Vq — So. Thus, the parameter /xi contains contributions coming 
from the scalar, as well as the vector part of the potential. The parameter /x corresponds 
to the rest mass of the Dirac particle. Substituting Eqs. (22), (23) and (24) into Eq. (7), 
we get 

d 2 k(k — 1) / e~ 2ar e~ 2ar \ 

x C"o " / VL,! + {V ° ~ go) e -^-l " Enk ) ) Gn ^ r) = ° ' (25) 
The following approximation is used for 1/r 2 term [31, 52] 



We compute 1/r 2 expanding into series. This will provide a physical result for a = 0. 
Defining a new variable s = e~ 2ar , we have 



1 - s d 1 
+ 



s(l-s)ds [s(l-s)] 2 

- 2r] 2 [fi (S - no) - E nK V + E 2 nk + 2a 2 K(n - 1) + /Ui/i ]s 

- Tj 2 [/ig - E 2 fc + 2£ nK \/ - 2SoA*o + m(m + 2S - 2fi ) 

- V 2 + S 2 ]s 2 })G nK (s) = 0, 



(27) 



where rf = I /4a 2 . Comparing Eq. (27) with Eq. (8), we get the parameter set given in 
Section II: 



01-1 = 1, 

a 2 = l, 

"3 = 1, 

0:4 = , 

"8 = 6, 



1 

4 ' 



-6 = V 2 [^o - E 2 n k + 2E nK V - 2S fi + + 2^0 - 2/i ) - V 2 + S 2 

6 = -2r] 2 [fi {S - /i ) - £™H + E 2 nk + 2a 2 /t(K - 1) + n m ] 

-& = v 2 (E 2 nk -i4) 

«7 = -6 

"9 = 6-6 + 6 + 1 



(28) 



a 10 = 1 + 2^, a u =2 + 2( ^ - 6 + 6 + \ + VTs ) 

«i2 = V6, «i3 = -1 - ( ^6 - 6 + 6 + I + v^s ) 

We can easily obtain the energy eigenvalue equation of the generalized Hulthen potential 
for any k value from Eq. (15) 



2^/io - E 2 nK = ^{2E nK V + fi 1 (fi 1 + 2S -2 f i )-V 2 + S 2 -2 f i S )-N, 



where 



(29) 



N = l -(2n + 1) + J V 2 \Vi(vi + 2S ) - V 2 + S 2 + Aa 2 K(K - 1)] + ^ . 



(30) 



The numerical results for different quantum numbers (n, /t) are listed in Table I. We list 
the eigenvalues E nK for constant mass = 0) and for two different values of ji\ to see the 



effect of the spatially dependent mass. The results given for constant mass are compared 
with results reported in the literature. For the constant mass case, we see good agreement 
with the results given in Ref. [42], where the relativistic energy W is stated as W = E + m. 
The lower spinor component can be obtained from Eq. (19) 



G nK (s) = S^(l - s y/^p(2e,2S) {l _ 2s) j 



(31) 



where e T 



rj 2 (pt - El) , and 5 = MM^i + ^ ) - K, 2 + S 2 + A(x 2 k{k - 1)] + \ . 



From the last equation, and using Eq. (5), the upper spinor component can be expressed 



as 



F nK (s) 



{ 



S .1 K 



H(s) - E nK + S(s) 
\{n + 2e nK + 25 + l)P 1 2+ 26 -' 1+M )(l - 2s)} . 



(32) 



We consider the case of constant mass to discuss the compatibility of our results. Setting 
Hi = in Eq. (24) while keeping in mind that V+(r) — C — const, for pseudospin symmetry, 
and following the same procedure (s = e~ 2ar ), we obtain the energy eigenvalue equation from 
Eq. (7) 



(^0 + E nK )(n - E nK + C) 



1 

+ rfV 



i(2n + l) + ^«(«-l) + ^ 
A*o — E nK + C 



i( 2n + l) + + i 



(33) 



where V = Vq + Sq . The result in the case of constant mass is Eq. (47) of Ref. [43] for 
So — > under the exact pseudospin symmetry. 

We now briefly give the bound state solutions for the constant mass case under spin 
symmetry. Setting V_(r) — C — const., and jii = 0, we get the bound state solutions under 
the exact spin symmetry for the case of constant mass from Eq. (6) as (s = e~ 2ar ) 



(fj, - E nK )(no + E nK — C) — 



An 2 
+ rfV 



\^n + l) + Jn{K + l) + - A 



C 



\ (2n + 1) + ^k(k + 1) + I 



where V' = S -V . 



(34) 
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V. CONCLUSION 



We approximately solved the Dirac equation, with position-dependent mass, for the gen- 
eralized Hulthen potential with arbitrary spin-orbit quantum number. We found the eigen- 
value equation, and corresponding two-component spinors in terms of Jacobi polynomials by 
using the parametric generalization of the NU-method within the framework of an approx- 
imation to the k(k — l)/r 2 term. We compared the numerical results with those obtained 
in the literature and given in Table I. We showed results for the case of constant mass, and 
summarized the results for two different position-dependent mass values, obtained when 
Hi = 0.005 and /ii = 0.0001. We also obtained the energy eigenvalue equation for the con- 
stant mass case with spin and pseudospin symmetries, separately. These analytical results 
are in agreement with results reported in the literature. 
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TABLE I: Energy eigenvalues for different values of n and k (/xq = Vq = 1) 



Mi 



a 



n k 



state 



E nK <0 



Ref.[47] 



0.025 
0.01 
0.05 
0.1 



1 -1 



Is 



1/2 



0.998068 
0.999691 
0.992258 
0.968772 



-1.993900 
-1.999000 
-1.975800 
-1.905700 



0.005 



0.025 



1 -1 
1 -2 
1 -3 
1 -4 



Is 



1/2 
lP3/2 



Id, 



5/2 



1/7/2 



0.998068 
0.995895 
0.993424 



0.05 



1 -1 

1 -2 
1 -3 
1 -4 



1*1/2 
lP3/2 
1^5/2 
1/7/2 



0.994106 
0.987615 
0.979641 
0.969946 



0.1 



1 -1 
1 -2 
1 -3 
1 -4 



1*1/2 
1P3/2 
1^5/2 
1/7/2 



0.970550 
0.945421 
0.912747 
0.871793 



0.0001 



0.025 



1 -1 
1 -2 
1 -3 
1 -4 



1*1/2 
1P3/2 
1^5/2 
1/7/2 



0.998103 
0.996592 
0.994653 
0.992281 



0.05 



1 -1 
1 -2 
1 -3 
1 -4 



1*1/2 
1P3/2 
l<4/2 
1/7/2 



0.992293 
0.986235 
0.978413 
0.968795 



0.1 



1 -1 
1 -2 
1 -3 
1 -4 



1*1/2 
1P3/2 



Id 



5/2 



1/7/2 



0.968807 
0.944015 
0.911464 
0.870558 
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